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Abstract 

In this paper, we study the performance of non-fading and Rayleigh fading ad hoc networks. 
We hrst characterize the distribution of the signal-to-interference-plus-noise ratio (SINR) through the 
Laplace transform of the inverted SINR for non-fading channels. Since most communication systems 
are interference-limited, we also consider the case of negligible noise power, and derive the upper and 
lower bounds for the signal-to-interference ratio (SIR) distribution under both non-fading and fading 
cases. These bounds are of closed forms and thus more convenient for theoretical analysis. Based on 
these derivations, we obtain closed-form bounds for both the average Shannon and outage rates. We also 
leverage the above results to study partial fading ad-hoc systems. These results are useful for investigating 
and comparing hfth-generation communication systems, for example massive multi-antenna and small¬ 
cell networks as in our illustrative example. 
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I. Introduction 

Due to frequency reuse, a wireless communication system is affected by not only the thermal noise but 
also interference from other transmitters in the same frequency band. The thermal noise can be modeled 
as additive white Gaussian noise (AWGN). On the other hand, the network interference characteristic is 
often elusive as it depends on many aspects, for example the spatial structure of the overall network or 
the propagation characteristics. For ease of analysis, interference can be approximated by AWGN with 
some given power IH. However, such simple model does not capture many physical parameters affecting 
the interference. More realistic models, using hexagonal or square lattice, are often difficult to analyze 
and thus mostly rely on simulations for investigation. System modeling based on stochastic geometry can 
account for all practical parameters while yielding explicit performance expressions, hence is promising. 

Stochastic geometry has been applied to investigate communication networks as early as 1970s l^. 
It has been successfully integrated and adapted to cellular, ultra-wideband, cognitive radio, relay, and 
especially ad hoc networks Q-Q. Since wireless systems is mostly interference-limited, a majority of 
the works have focused on characterizing the interference and the signal-to-interference ratio (SIR) at each 
receiver/user equipment (UE). Please refer to lITOl - lfT^ for detailed summaries. Under most networks, 
the interference and/or SIR distribution were represented only through Laplace transforms, which are 
often complicated. Useful bounds of the SIR distribution for several system models have been proposed 
in IIT3I and lIT^ . In the seminal work lITSl . the fading ad hoc network has been investigated under the 
assumption that the desired signal power at a typical receiver/UE is exponentially distributed. Based on 
the Laplace transform and probability generating functionals of Poisson point processes (PPPs), Andrews 
et. al. have shown that the interference and SIR distribution of a wide range of networks can be greatly 
simplified. Such approach has enabled simpler derivations of the outage probability and achievable rate, 
and thus motivated more investigations llT6ll - ll^ . 

Stochastic geometry is particularly useful in modeling and studying heterogeneous networks, which 
consists of many coexisting communication systems and in which all participating transmitters and 
receivers are randomly located 11191 - 11211 . Coordination and cooperation between cells have also been 
considered, e.g., in ifTTl . A shortcoming of many existing works is that they assume a very large UE 
density compared to the transmitter/BS density, effectively considering severely interfering networks only. 
To address this shortcoming, several recent studies have introduced a general model which includes the 
UE density as a varying parameter for performance investigation lIThl . ifTSl . 

In this paper, we obtain new results for the signal-to-interference-plus-noise ratio (SINR) and SIR 
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cumulative distribution function (CDF) at a typical UE under non-fading and Rayleigh fading ad hoc 
networks, where the transmitting nodes are distributed according to PPPs. Particularly for non-fading 
ad hoc network, the Laplace transform of the inverted SINR are given, with which the distribution can 
be numerically calculated. We also derive upper- and lower bounds for the SIR CDF, which possess 
closed-form expressions and are thus more useful for subsequent analyses. For Rayleigh fading ad hoc 
network, we obtain closed-form bounds for the SIR CDF. Based on these results, the average Shannon 
and outage rate bounds are acquired under both cases. Similar results are obtained for partial fading 
ad-hoc systems. Since many communication networks can be modeled or reduced to these two ad hoc 
systems, our results are potentially useful for a wide range of problems. For example, we will show that 
our results can be applied to study and compare the performance of two important fifth-generation (5G) 
technologies, i.e., massive multiple-input multiple-output (M-MIMO) and small-cell densification. More 
information can be found in our subsequent paper |[22l . 

The rest of the paper is organized as follows. Section HIl introduces the ad hoc network and its system 
model based on stochastic geometry. In Section |III1 we derive the Laplace transform of the inverted 
SINR, and the bounds of the SIR distribution and rates for non-fading ad hoc networks. Similar results 
for Rayleigh fading and partial fading ad hoc systems are obtained in Sections |IV] and |Vl respectively. 
Section |Vl] presents one application of our results to study M-MIMO performance and to compare it with 
that of small-cell densification. Finally, Section IVIII concludes the paper. 

Notations: Scalars and vectors/matrices are denoted by lower-case and bold-face lower-case/upper-case 
letters, respectively. The conjugate, transpose, and conjugate transpose operators are denoted by (•)*, (-j^, 
and i , respectively. [A]ij stands for the (i, j)th element of the matrix A. Ex[-] denotes the statistical 
expectation over a random variable X. Tr(-), || • ||F^ and det(-) represent the trace, Frobenius norm, and 
determinant of a matrix, respectively. Im denotes an M-hy-M identity matrix. and denote 

the space of x-hy-y real and complex matrices, respectively. Finally, the circularly symmetric complex 
Gaussian distribution with mean y and variance is represented by 

II. System Model eor Ad-Hoc Networks 

In this paper, we consider ad hoc wireless networks with an infinite number of base stations (BSs) 
randomly distributed over an infinite area. The locations of the BSs are modeled by a homogeneous 
PPP <1> C with density A. The single-antenna UEs are assumed to be distributed according to 
some independent stationary point process, and each UE is associated with the nearest transmit node. 
Furthermore, the UE density is large enough such that the BSs all transmit at any given time. The resulting 
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system model is the so-called ad hoc network in the literature ifTOll - llTSl . Consider an arbitrary BS 6. 
Under the power-law path loss model, the hase-hand received signal for the UE in cell b is given as 

OO 

Vb = hbfid^l^^'^Xb + ^ hb,kdl^^^‘^Xk + rife, (1) 

k^b 

where db^k and /ife^ are the distance and the independent and identically distributed (i.i.d.) fading 
coefficient between BS k and the UE associated with BS 6; x^ is the transmit signal of BS k, where 
E [|xfcp] = Pt is the transmit power; Ub is the additive white Gaussian noise (AWGN), Ub ~ CA/"(0, ci^); 
and /i is the path loss exponent. 

Since the system model is homogeneous, all analysis for a typical UE is applicable to any other user 
due to the Slivnyak’s theorem |[T0l - llT2l . Eor ease of notation, hereafter we consider the UE associated 
with BS 1 and drop the index 1 for brevity. The SINR for a typical UE can thus be given as 

\hifd^^ 


SINR = 


( 2 ) 


EZ2\hkf d->^ + 6' 

where 6 = /Pt is the normalized noise power. Eig. [T] shows a channel realization of an ad hoc network 
with BS density A = 0.05 km“^ and service area assumed to be 20 x 20 km^. 


III. Performance of Non-Eading Ad-Hoc Networks 

In this section, we consider a non-fading ad hoc network where |/ifcp = 1, \/k, in the system model 
described in Section |II1 Such system is also referred to as the base-line model |[T3l . We first obtain 
the Eaplace transform of the inverted SINR. However, the resulting expression is complicated and thus 
difficult to use for theoretical analysis. We then derive the upper and lower bounds for the SIR distribution, 
which are in closed forms. Based on these results, we investigate and obtain the bounds for both the 
average Shannon and the outage rates of non-fading ad hoc networks. 

Consider the random variable 


Qn-fd — 


d'l 


-n 


(3) 


Er=2 4"+^ 

It represents the SINR at a typical UE under non-fading ad hoc networks with noise power 5. However, 
existing works have assumed that the distance d\ is a constant and only the set {dA :}^2 follows a PPP 
ma, m- in contrast, we consider the case when follows the PPP with density A and dk> di, 

\/k > 2. This is a typical regime where the UE is associated with the nearest transmitter. The distance from 
the UE to the associated transmitter, di, is therefore smaller than those to other transmitters, {dfc}^ 2 - 
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Fig. 1; A channel realization of an ad hoc network. 
Here, and straight line denote the BSs, UE, 

and boundaries between the Voronoi cells. Furthermore, 

and “-” denote the channels between UE and 

its associated and interfering BSs, respectively. For ease 
of presentation, only one UE in one cell is shown here. 


We first derive the Laplace transform of (5n-fd^ based on which the distribution of Qn-fd can be calculated 
numerically. We furthermore propose CDF upper and lower bounds of Qn-fd^ which is much more useful 
for theoretical analysis than the Laplace transform. 

Theorem 1: The Laplace transform of the random variable is given as 


^l/Qn-a (^) - 


exp — sy^^‘^5 — nXy 

2 — ^1 — —, s 


dy, 


(4) 


Proof: Please refer to Appendix |3 ■ 

Theorem [U provides a semi-analytical approach to calculate the distribution of and consequently 
of Qn-fd- However, such a characterization requires intensive numerical computations, but does not provide 
any insightful observation. We thus seek to find bounds for the CDF of Qn-fd which is more useful for 

d-P 

theoretical analyses. Here we assume that the noise power 5 is small and approximate Qn-fd ~ - 

^fc = 2 

In M-MIMO systems, this approximation is well justified due to effect of M-MIMO beamforming which 
compels the noise power to much smaller than the inter-BS interference. The following lemma establishes 


September 15, 2015 


DRAFT 










6 


upper and lower bounds of the CDF of Qn-fd given that (^ = 0. 

Lemma 1: Given that (5 = 0 and q > 1, the CDF FQ^f^{q) = PrjQn-fd < q} is upper- and lower- 
bounded as F^lJq) < FQ„,,(g) < F^lJq), where 


^qL(^) = i- 


1 


1 


/3n-fdg^/^ ’ 


(5) 

( 6 ) 


and 


_ (p -h 2)2//"+^ 
fj,{p — 2 ) 2 //^ 


(7) 


Proof: Please refer to Appendix |B] ■ 

It is difficult to investigate the regime (7 < 1. In the following, we provide an asymptotic upper bound 
for Qn_fd as g —)• 0 or g — 1 . 

Lemma 2: Given that (5 = 0 and q <1 but g —0 or 1, the CDF FQ^_i^{q) is upper-bounded as 

1 


^jn-tc 


( 9 ) = 1 - 


1 + /3n-fdg'^/^ ’ 


( 8 ) 


Proof: Please refer to Appendix O 

Based on Lemmas [T] and |2l we propose the following universal bounds for Qn-fd- 


FqIM) 


0 , 

1 

g2 /fi ’ 


0<q<l, 

q > 


(9a) 

(9b) 


7^Qi(7) = l- 


1 + /3n-fdg^/^ 


( 10 ) 


Analytically, (fTOl ) has been proved only when: ( 7 >l;( 7 ^ 0 ;or( 7 <l but q 1 in Lemmas [T] and |2] 
However, intensive simulations show that it holds for all values of q. In Fig. |2l we show the simulated 
SINR CDF and the bounds given in ®, (|9l), and (fTOl) . for an ad hoc network with /j, = 3.7, S = —20 
dB, and 20000 channel realizations. We see that the bounds are close to the exact distribution. With a 
path-loss exponent following that of LTE, i.e., p € [3.7,4] (241, the difference between ® and the looser 
(but universal) (fTOl) is negligible. However, as p increases, we observe that (fTOl) deviates farther from the 
exact CDF than ®. Another observation is that the tail of the exact SINR CDF seems to follow a mixed 
scaling law ~ x““exp(— 6 x). This observation might help subsequent works for developing tighter CDF 
bounds. 

The following Corollary [T] gives the expectation of averaged over the PPP <h. This result might 
be useful for comparison between different systems. The proof can be deduced from Appendix |Cl A 
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Fig. 2: Simulated SINR CDF and its bounds for non¬ 
fading ad hoc network with fj, = 3.7 and S = —20 dB. 



Fig. 3: Average of with respect to the path loss 
exponent. 


confirmation for Corollary [T] is illustrated in Fig. |3l where we plot the inverted SINR versus the path-loss 
exponent with (5 = — 20 dB 0. 

Corollary 1: For non-fading channel, the expectation of the inverted SIR is given as [Q~ = 0] = 

2 

fi—2 ■ 
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A. Average Achievable Rate 


In this subsection, we obtain several results for the average Shannon rate of a typical user, defined as 


-Rn-fd(a) = IE$ [log 2 (1 + aQn-fd)] , 


( 11 ) 


where Qn-fd is given in ([3]) and a is a constant representing, e.g., the beamforming gain of the system. 
Essentially, we have assumed that single-user receiver is employed, i.e., the interference is treated as 
noise. We first note the following fact 

Fact 1: For a positive random variable X, the expectation of X can be expressed as E [X] = Pr {X > t} dt. 
Based on Fact [T] the average rate can be written as 


poo 

-Rn-fd(a) = / Pr {log 2 (1 + aQn-fd) > t} 
Jo 


dt 


f 


Pr S Qn-fd ^ 


2 * - 1 


a 


dt. 


( 12 ) 


-1 

fd 


Therefore, an exact i?n-fd(a) can be calculated numerically by using the Faplace transform of Q~ 
in Theorem [T] Such approach, however, does not give us any insight on the system performance. We 
therefore provide bounds of i?n-fd(a) based on ® and (ITOl) in the next lemma. 

Lemma 3: The average Shannon rate i?n-fd(a) is upper- and lower-bounded by and a), 

respectively, where 


-^n-fd(rr) — log 2 (l -|- a) -|- 

= [ 

Jo 




(p - 2) log 2 
o?l^dt 


P a 


1 - 2 //. _J^ 


’ o / ’ 

/i — 2 / 


0 a^/t*-h/3n-fd (2* - 1) 


2/ti- 


(13) 

(14) 


A dx 

Ja x‘’+t 

Proof: Please refer to Appendix iDl 


where p(o, h) = 


B. Outage Rate 

Outage rate is another performance metric which for our system can be defined as 

ORn-fd(a, r]) = Pr |Qn-fd > log 2 (1 + ^), (15) 

where a is a given beamforming gain and rj is the target SINR corresponding to a target rate of log 2 (l-|-? 7 ). 
Similar to the Shannon throughput, the exact outage rate can be numerically determined based on the 
Faplace transform ([D). Again, we derive the bounds for the outage rate in the following lemma. 
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Fig. 4; Average Shannon rate and its bounds ( fTSI l. (fTTb 
for non-fading ad hoc network with fi = 3.7. 


Lemma 4: The outage rate ORn-fd(a,r?) with a > 0 is upper- and lower-hounded hy ORj^_f£i(Q!, 77 ) and 
a,ri), respectively, where 


OR^l{a,v) 


'log 2 {l + r]), 0<r]<a, 

log 2 (l+ r?) ,r]>a. 


(16a) 

(16h) 


ORn%{a,rj) 


log2 (1 + T?) 


(17) 


where /3n-fd is defined in Lemma [T] 

Proof: The proof is obtained hy substituting ([5]) and (fTOl) into (fTSl) . It is thus omitted for brevity. ■ 
The outage capacity of a system can be defined as the maximum of ORn-fd(a, 7 ) over rj, i.e., OCn-fd(a) = 
max^>o ORn-fd(a, ??)■ However, it is not easy to derive OCn_fd(Q:) based on Lemma ID 

In Figs. |4] and [5J we illustrate the average Shannon rate and the outage rate over a and rj, respectively. 
For comparison, we also plot their analytical upper and lower bounds given in Lemmas [3] and |4l It is seen 
that the upper bound is tighter than the lower bound for the average rate. The observation is reversed for 
the outage rate, in which the lower bound is (slightly) tighter. One reason for the slackness of the lower 
bound is due to the overestimation of (ITOl) . Intensive simulations suggest that the tail of the exact SINR 
CDF follows ~ x““exp(— 6 x), which approaches 0 faster than the polynomial law given by the 

current CDF bounds (|9l) and (ITOl) . To obtain tighter CDF bounds, more stringent analyses are necessary. 
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Fig. 5; Outage rate and its bounds ( fOl l. (O for non¬ 
fading ad hoc network with respect to the target SINR. 
Here fj. = 3.7 and a = 10 dB. 


Remark 1: Assume that the noise power is ignored, i.e., 5 = 0. From TheoremlH we observe that the 
Laplace transform of does not depend on the density A, and it can he expressed as 

2 - e"^ 

As a result, the CDF and rate of non-fading ad hoc networks are independent of the density. Such 
observation is readily demonstrated in Lemmas [T] to |4l At first glance, it seems to contradict with the 
results in ifTSl . ll23l . The reason for the difference is that Weber et. al. only assumed the interfering BS 
locations to distribute according to a PPP, but not the associated BS. The distance between the UE and 
associated BS is a constant in ifTSl . However in this study, we assume that all BS locations are 
distributed according to a PPP. 

We also note that the independence of the CDF and rate to the density can be explained geometrically. 
Having a higher density resembles the effect of zooming in a certain region and adding more transmitting 
nodes. Since the considered PPP is homogeneous, the resulting PPP is stochastically equivalent to the 
original one. The relative ratios between the distances are therefore unchanged. It explains why the CDF 
and the rate are unaffected by the BS density. 
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Fig. 6: Simulated SINR CDF and its bounds (1211 1. (l22l) 
for fading ad hoc network with /j, = 3.7 and h = —15 
dB. 


IV. Performance of Fading Ad-Hoc Networks 


In this section, we consider a slightly different prohlem than Section [Till Particularly, we investigate 
the quantity 


Qfd = 


\hi\U 


2 


(19) 


where <I>, {dk}’^=i, and 6 have been defined in Section ini In the literature, Qa is the SINR of 

an ad hoc network under fading channels, and also attracted much attention. The exact outage prohahility 
of Qfd is given in ifTSl for the general and various special cases. The crucial assumption in IfTSll is the 
exponential distrihution of \hif, which allows a simple representation of the SINR distribution. 

We are interested in the Rayleigh fading channel, i.e., is exponentially distributed with variance 
1, of which the SINR CDF is provided in ifTSl Theorem 2]. Assuming that the noise power is h = 0, a 
much simpler expression can be obtained as lITSl 


FQiA<l) 


Pr {Qfd < q} = 1 - 


1 

1 + ’ 


( 20 ) 


where the function p(a, b) is defined in Lemma [3] The expression (l20l) requires a numerical integrafion 
and is thus not useful for performance analysis. We will develop the lower and upper bounds for (l20b . 
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Lemma 5: Assume (5 = 0. The lower and upper bounds for the CDF of Qfd defined in ([T^ are given 


as 


^q"(9) = i- 


1 29=^" p 


2 

Q 


where h ~ CJ\f{Q, 1); and 


/3fd = 


1 + 


2tt 


=r(1+- 


2 

1 - 

F 


( 21 ) 

( 22 ) 

(23) 


H sin 

Proof: Please refer to Appendix 10 ■ 

Fig. [0 depicts the SINR CDF and the bounds given in Lemma[5]for ^ = 3.7 and 6 = —15 dB. Unlike 
the non-fading network (see Fig. |2), the bounds for the fading case is much stricter, especially the upper 
bound. For path-loss exponents under typical LTE regimes, i.e., ^ = 3.7 to 4 ll2^ . the upper bound is 
very close to the exact CDF, and gradually becomes looser as p increases. It is expected that the gaps 
between the exact rates and bounds given in Lemma will be smaller than those under the non-fading 
case. 

In the following, we derive the average achievable Shannon and outage rate bounds based on the CDF 
upper bound. The proofs are similar to those in Sections ITTT-AI and ITTT-RI and thus omitted for brevity. 

Lemma 6: The average Shannon rate i?fd(a) and the outage rate ORfd(Q:, ??) are bounded above, 
respectively, by 


R\i(a) = 


I 


OR\^{a,T]) = 


0 + 

log2 (l + v) 


(24) 

(25) 


The rate upper bounds can also be derived, however, they possess complicated integral forms with 
little practical use. The reason is due to (l2TI) which requires numerical calculations. We therefore do not 
present the rate lower bounds here. 

In Figs. |2] and [H we illustrate the average Shannon and outage rate over a and rj, respectively, for 
Rayleigh fading ad hoc networks. For comparison, we also plot their analytical lower bounds given in 
Lemma [0 and the analytical values calculated based on the noise-free CDL approximation (l20l) . As 
expected, the bounds for the fading channel are much tighter than the non-fading counterpart given in 
Section 1111] 
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Fig. 7; Average Shannon rate and its bound (l24l l for 
Rayleigh fading ad hoc network with /i = 3.7 and S = 
— 15 dB. The analytical values are calculated based on 
the noise-free approximation (l20l l. 


Similar to the non-fading case, here the CDF and rates are also independent of the density. The reason 


is the same as in Remark [T] Finally, we note that given <5 = 0, the average of the inverted SIR E [Q 




goes to infinity since E/^ 


\h\ ^ is unbounded, where the random variable h ~ CA/’(0,1). 


V. Performance of Partial Fading Ad-Hoc Networks 

In this section, we will investigate partial fading ad-hoc networks, under which only the interfering 
signals experience Rayleigh fading channels while the desired signal does not. The SINR of a typical 
user is expressed as 

^P-fd = ,! ,2 1-tl I e ’ 

where <1>, {dk}’^=i-> and <5 have been defined in Secfion|ni In the literature, Qp-fd has also been 

studied in ifTSll . ll23l . However, as discussed in Section Hill the distance di is assumed to be a constant 
and only the set {dfc }^2 follows a PPP ifTSll . 12^ . Here, we consider the case when follows 

the PPP <1> with density A and > di,\/k > 2. 

Similar to Section |III1 we will derive the Laplace transform of and the distribution bounds for 
Qp-fd assuming that <5 = 0. The results are presented in the following lemmas. Note that no upper bound 
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Fig. 8; Outage rate and its bound ( l25l l for Rayleigh 
fading ad hoc network with respect to the target SINR. 
Here = 3.7, a = 10 dB, and 5 = -15 dB. We also 
show analytical values calculated from the noise-free 
approximation (l20l) . 


is available due to the difficulty of the derivation. 


Lemma 7: The Laplace transform of is given as 


'^i/Qp-fd (s) = 27rA xexp <* - 


— vrA 


+ 2 


1 - 


1 


dv 


dx 


1 + sx^^v~f^ 

Proof: Please refer to Appendix |F] 

Lemma 8: Assume (5 = 0. The lower bound of the CDF of Qp_fd defined in 


(27) 


IS given as 


pm 


Qp-fd(7) - 1 - 



(28) 


Proof: The proof follows from (ISTI) in Appendix |El ■ 

It is difficult to derive a CDF upper bound of partial fading ad-hoc networks. In the following, we will 
present an interesting relationship between non-fading and partial fading networks which leads to such 
bound. Note that given a PPP configuration, ^ summation of an infinite number 

of independent random variables. Without loss of generality, we can rank {dfcjfcg# as di < d 2 < ■■■■ Due 
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Fig. 9; Simulated SINR CDFs of Qp-fd and Qn-fd and 
their bounds with fj, = 3.7 and <5 = —20 dB. 


to the Kolmogorov’s strong law of large number Il25l Theorem 2.3.10], we have 


N r N , N 

|2 








k=2 


k=2 




(29) 


k=2 


Since 


k=2 k=2 

with high probability when |U > 2. Therefore, we can approximate 


(30) 


E E 4'‘' 

fce$\{i} fce#\{i} 

and Qp-fd ~ Qn-fd- With this approximation, we can utilize all results, not only restrict to the upper bound, 
from Section |III] to study partial fading ad-hoc networks. 

In Fig. |9j we show the simulated SINR CDFs of non-fading and partial fading ad-hoc networks as 
well as their bounds. Here, the number of channel realizations is 20000. It is observed that the CDFs of 
Qp-fd and Qn-fd are closely matched. Furthermore, even through simpler, the lower bound (|28] ) is much 
looser than ® Together, those observations well justify the approximation Qp_fd ~ Qn-fd and the use of 
the results in Section |III] to study Qp-fd- 


September 15, 2015 


DRAFT 



























16 


VI. Case Study: Massive MIMO and Small-Cell Densification in 5G 


Towards the development of the fifth generation (5G) communication systems, there emerge many 
promising technologies. It is therefore necessary to compare their performance under different scenarios, 
developing insights to deploy the right technology in the right scenarios. Note that among the technologies 
are M-MIMO and small-cell densification, both promising higher achievable throughput with either 
concentrating massive number of antennas in one location and serving the users with narrow beams, 
or distributing the large number of antennas close to the users’ proximity. 

The results obtained in Sections |III] and |IV] are useful to obtain performance insights for 5G systems, 
since many of them can be modeled or reduced to non-fading and Rayleigh fading ad-hoc networks. 
Furthermore, it is possible to compare the performance of such systems using the bounds given in 
Sections |III] and |IVl due to their closed-form expressions. 

Here, we give an example of M-MIMO networks. Consider the system model in Section |III] but with 
M antennas at each BS. The base-band received signal for the UE in cell b is given as 


CX> 

Vb = ^ + nb, (32) 

b'^b 

where h^, ^ is the channel between BS b' and the associated UE of BS 6; and Xb' is the transmit signal 
from BS b'. Here, hb'^b = ib',b - d^i^b ^ large scale fading including path loss with /r 

denoting the path loss exponent; hb'^b follows a statistical distribution hb'^b ~ CA((0,1); and Ub is the 
additive white Gaussian noise (AWGN) distributed as CA((0,cr^). 

Assume that each BS has a total power constraint of MPt and therefore E|a;b/p = MPt- Furthermore, 
assume that each BS employs the (optimal) conjugate beamforming, i.e., Xb = WbSb = where 

Wb and Sb are the beamformer and desired signal of BS b, respectively. We can express the SINR of the 
UE associated to BS 6 as 


SINRf, 


db,i: 


h. 


b,b 


/M 


-H 






u ^b',b^b',b'^b',b'^b',b 


M 


hb',b' 


+ 


cr„ 


M^Pt 


M—>-oo 
-^ 






(33) 


where hb' = 





are i.i.d. random variables each distributing as CA^(0,1). 
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It can be observed that (l3^ has a similar form to (|2^ . We thus can apply the results in Sections ITm and 
Ivlto investigate the performance of M-MIMO networks modeled as in (l3^ . Note that small-cell systems 
can be modeled as Rayleigh fading ad-hoc networks. Therefore, we can further compare the performance 
of M-MIMO and small cell networks due to the closed-form expressions of the performance bounds in 
Section JII] and |IVl Please refer to our subsequent study Il22l for more details. 

VII. Conclusions 

In this paper, we have investigated the non-fading and Rayleigh fading ad hoc networks. In particular for 
the non-fading case, we have derived the Laplace transform of the inverted SINR, which is unfortunately 
complicated and of limited use. Therefore, we have also proposed upper and lower bounds for the SIR 
distribution and given the closed-form expressions. The results have then been used to obtain the bounds 
for both the average Shannon and outage rates. Similarly for Rayleigh fading and partial fading ad hoc 
networks, we have derived the bounds for the SIR distribution, average Shannon rate, and outage rate 
of the system. Since many communication networks can be modeled or reduced to these two ad hoc 
networks, our results are potentially useful for a wide range of problems. As an example, we have 
shown that our results can be applied to study and compare the performance of two 5G communication 
technologies, i.e., massive MIMO and small-cell densification. 
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Appendix A 
Proof of Theorem [H 

We first derive the Laplace transform of the random variable given di = x as follows. 


{s\di = x) = E^ [e-^Q"-“] 

= ]E$ 


CXD 


exp I —sx'^ ^ — sx^6 

k=2 / 

OO 




(®) —sx^5 

= e exp 




n 

k=2 

{-n 


I — e 


— finrJ^r 


2'K\rdr 


0 ^_,^.5exp<; - t:\x\ 1-e-^) 


J X 

(J ^-sx-& I _ _ g-.) 


(^) —sxf^S 


= e 


exp 


-ttXsx^ [ 

Jo 

|—7rAx^(l — e~^) — 7rAs^'^^x^7(l — 2//i, s)| 


(34) 


where (a), (6), (c), (d) is due to the Camphell’s theorem IIT^ . the integration hy part, the change of 
variable v = r~^, and ll26l (3.351.1)], respectively. Here 7 (-, •) is the lower incomplete gamma function, 
defined as 'j{a,b) = t°‘~^e~*dt. 

Now note that di is the distance from an arbitrary origin to the nearest point of the PPP <h with density 
A. Therefore, the PDF of di is given as ifTSl iflTl 


fo, {x) = 27r Axe . 


(35) 
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The Laplace transform of is thus expressed as 

('S) 

roo 

^Qn-m = x) fD^{x)dx 


IQ 


= 27tX J X exp — sx'^S 


— ttAx^ 

= ttA J exp I — sy^^‘^6 

- vrAy 


2-e-" + s2/M7(l--,s 


2 — e ^ ( 1-, s 


where (e) is due to the change of variable y = x‘^. 
This concludes the proof of Theorem [T] 


dx 


dy, 


(36) 


Appendix B 
Proof of Lemma [T] 

A. Lower Bound 

Considering q>l and the second smallest distance d 2 only, we have 

FQn-tdiQldi = x) > Pr 1^2 < xq^^^\di = x| 


(a) 


rxq 


1/fU. 


J X 


(37) 


where (o) is due to |[27l (28)]. The PDF of di follows (1351) . which thus leads to 

roo 

= x)fDdx)dx 

Jo 

= (l - 2'KXxe-^^^"dx = 1 - (38) 

B. Upper Bound 

Assuming that q > I and ^ > 1 is an arbitrary parameter, we now apply the dominant-interferer 
concept lIT^ . The idea is to partition the set of interferers into dominating and non-dominating points, in 
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which a point /c (> 2) is called dominating if it alone can cause outage at the UE. The sets of dominating 
and non-dominating interferers, denoted as <Sdom and ^Sn-dom^ respectively, are defined as follows 


>5dom = : k> 2, dk < , 

*5n-dom — /c . k>2,dk> . (39) 

For ease of representation, we also define Sn_dom = '^kes„d m ^idom = X]fce<Sdom 

that 


Qn-fd — 


d- 




^fe65d„m ^ + X]fce<S„.do„ 


-M 


< 


d-^ 


Sfc 65 d„„ ^ + Y.k&Sn-dom ^k ^ 


-u ~ Qn- 


fd- 


We thus have 


(40) 


Fq..M < Pr {4-fd <q}= F^^Jq). 


Given that di = x, we have 


Pr 


|Qn-fd < 9|c^i = a^} 


g 


— Pr <; <; Sdom > - J U | Sn-dom > 

U \ Edom ^ ~ 5 ^n-dom ^ ~ j ^dom T ^n-dom ^ ~ 

q q q 


di = X 


(“) n 

= Pr 


{b) 


1 


-'dom 


q 


1 


> - ^ U Sn-dom > - 


q 


= Pr <! Sdom > - 

q 


d\ — X ^ T Pr Sn-dom ^ ^ 


di = X 
1 


di = X 


- Pr Edom > 


dl — X > Pr s Sfi-dom ^ 


= Pr Edom > 


di = X ^ + 1 - Pr Edom > 


di = X 
1 


q 


di = X 


/e > ^ 

di = X 

1 q 

J 


(41) 


(42) 


where (a) is due to the fact that the event {Edom < l/ 9 |di = x} is null; (b) is because of the independent 
of 5dom and iSn-dom- From (iTTl) . we obtain 

1 


Pr { Edom > 


di = x| = Pr |d 2 < x(^( 7 )^'^^|di = x| 
= 1 _ 


(43) 


September 15, 2015 


DRAFT 














21 


Furthermore, we have 


J 1 

) 

(c) 


S Efi-dom ^ ~ 

1 Q 

di = X 

Crf 

VI 

^n-domMl — ^ 


(44) 

where (c) is due to the Markov inequality Pr{X > a} < Now hy applying Camphell’s theorem, 

we further obtain 


E# [Sn-dom|dl = X] = 


E < 


k 


f 


V ^27rXvdv = 


keS„ 

-dom 


Substituting (|45] ) and (1431 ) to (1441) . the upper bound given that di = x is expressed as follows 

FZ.m = X) = 1 - 

fi — 2 ’ 


(45) 


(46) 


and consequently we have 


POO 

= x)fD^ix)dx 

Jo 

= 2vrA^°° (l - xe-^^^"dx 



+ 


{2xx)^e/^-v/>^ 

Jo 


fi — 2 


1 2^2/At-lg2/M I 

(^g)2/M // - 2 (^q)^//" 

(^ - 2)^ - 2 1 
{^l-2)e+^/^^ g2//.’ 


(47) 


where (e) is due to (l3^ and |[26l (3.321.6)]. The final step is to find the maximum of given 

that C > 1- The optimal is The result is ®. 

This completes the proof of Lemma [T] 
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Appendix C 
Proof of Lemma [2] 
Using the Markov inequality, we first have 

FQnA(l\di = x) = Pr > q-^\di = x} 

< qE^ [Qn_fdl*=^i = =^‘S> 


X‘ 


k>2 


= qx^ 


V ^2TT\vdv = — —x^. 

11 — 2 


(48) 


The CDF is thus hounded above as follows 

FQn-iAo) < [ = -^^q- (49) 

Jq 11-2 q-2 

Before proceeding, we note that -j^q does not give a strict hound for q > 1 or q < 1 hut q ^ 1, 
since 1 — (^)|) 2 ) 2 /a.+i < Now it is observed that, as q —)> 0, 


-q<l- 


1 


(50) 


11 - 2 " ' 1+ /3n-fdq^/^' 

Furthermore based on Lemma [H it is straightforward to show that 1 — is an upper bound 

when q —> 1. 

This concludes the proof of Lemma |2] 


Appendix D 
Proof of Lemma [3] 


We have 


r°° f 2^ — 1 
din-fd(a) = J Pr |g > —-— )■ dt 

/■logafl+o:) 


4o 


1 -Pr|Q < 

1 - Pr I g < 


Substituting (|9l) to (ISTI) . we have 


-Rn-fd < Rn-fd “ 


logafl+a) 


r'logjfl+a) poo 

dt I 

3 die 


2 * - 1 
a 

2 ‘ - 1 
a 


dt 

dt. 


a'^/f^dt 


- log2(l + a) + 


iogAi+«) (2* - 

dv 


1 — - ) log 2 4iog2(i+a) n 1 - 2 /m + 1 


^log,(l + a)+ , “'dr ,, 

^ (/x-2)log2^V V- 2 ;’ 


(51) 


(52) 
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where (a) is due to the change of variable v = (2* — l)^ The lower hounds can he obtained with 
similar steps by substituting (fTOl) into (ISTI) . and thus are omitted for brevity. 

This completes the proof of Lemma |3] 


Appendix E 
Proof of Lemma [5] 

A. Lower Bound 

Given hi and di, we first define 

>5fd.dom = |fc > 2 : \hkfdl^ > . 

The SINR CDF can be bounded below as 


(53) 


FQiMdi,hi) > FQ^^{q\di,hi), 


where 


FQliQ\di,hi) =Pr < 


l^il ‘^d^ \hk\‘^di^^ > q ^ 

fc€<5fd,dom 


—^ 1 — Pr {5fd,dom — 0} 


(b) 


(54) 


= 1 — exp ( — y Pr ||/ip|x| ^>q ^\hif d-^^ dxj , (55) 

where (a) is due to the fact that the event {<Sfd,dom ^ 0} is the same as |X]fce<Sfddom \dk\‘^d^^ > q~^\hi\^d^^'^', 
and (6) is due to the void probability of a Poisson process. Using the polar coordinates and the fact that 
|/ip is exponentially distributed, we have 


Pr||/i|2|xr^ > 


dx 




/•CO 

Jdi 


■IrrXr dr 


W 27rA ^ / 2 |/ii|2 


h \hi\yt^ \fr’ q 


(56) 
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where (c) is due to ll^ (3.381.9)]. The CDF lower hound thus can he obtained as 




exp 


2tt\ d\q^/^ ^ ( 2 |/iip 


= 1 - E;,, 


[/ 


exp 


IX |/ii|2/m q 

2g2/M ^f2 \hi\ 


X 27r\te~^^^ dt 


ttAC 


= 1 - E/i 


1 + 


2q^/i‘ 

TfTlTF 


.r I 2 Ihiil 

‘ I /i’ g 


(57) 


B. Upper Bound 
We note that 


f 

Jg 


du 


< 


du (d) 


27T 


2/m 


1 + tt^/2 Jq 1 _|_ ^m/2 




= r( i + - )r( 1-- 


where (d) is from |[2^ (3.241.2)]. As a consequence, we have 

^Qaid) — 1 - 2/fi TOO rfn 

1 

< 1 - 


l + /3fdy2//^ ■ 


This completes the proof of Lemma |5] 


(58) 


(59) 


Appendix F 
Proof of Lemma |7] 

We first derive the Laplace transform of the random variable Q~f^ given di = x as follows. 


El/Op.H (s|dl = x) = [e 


= E, 




exp ( —— sx^5 


k=2 








.k=2 


(a) 


= exp ( — s6x^ — 27rA 
where (o) is due to (HI (39)]. 


1 - 


1 


1 + sx^u ^ 


dv 


(60) 
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Again, since the PDF of di is given as in (1351) . the Laplace transform of is thus expressed as 


-^l/Qp-fd (s) 



(61) 


This concludes the proof of Lemma |7] 
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